Damping of bulk excitations over an elongated BEC - the role of radial modes 
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We report the measurement of Beliaev damping of bulk excitations in cigar shaped Bose Einstein 
condensates of atomic vapor. By using post selection, excitation line shapes of the total population 
are compared with those of the undamped excitations. We find that the damping depends on the 
initial excitation energy of the decaying quasi particle, as well as on the excitation momentum. We 
model the condensate as an infinite cylinder and calculate the damping rates of the different radial 
modes. The derived damping rates are in good agreement with the experimentally measured ones. 
The damping rates strongly depend on the destructive interference between pathways for damping, 
due to the quantum many-body nature of both excitation and damping products. 



Ever since Wigner and Weisskopf first calculated the 
damping rate of an atom coupled to the vacuum [1], it 
is known that coupling to a bath affects the decay of 
an otherwise stable quantum system. The structure of 
the bath plays a key role in determining the rate of the 
damping, and may even cause the damping to be sub- 
exponential or super-exponential Q]- Gaseous Bose Ein- 
stein condensates (BEC) are usually well isolated from 
their surroundings, leading to coherent evolution. Three 
body loss, which is the main cause for decoherence in the 
ground state, can be on the order of seconds However, 
damping of the elementary excitations over the conden- 
sate is much more rapid, but still accessible experimen- 
tally. The bath, which is coupled to the excitation, is 
in this case a quasi-continuum of unoccupied excitations. 
These excitations, known as Bogoliubov quasi particles, 
are approximate eigen states of the Bosonic many-body 
Hamiltonian. Coupling of an excited quasi particle (QP) 
to initially unoccupied QP modes gives rise to decoher- 
ence via the Beliaev damping mechanism, in which a 
QP decays into two new QPs while fulfilling energy and 
momentum conservation [J|. Another damping process, 
that takes place at a finite temperature, is called Lan- 
dau damping. This process involves the annihilation of 
a thermally excited QP together with the damped one, 
and the creation of a more energetic QP, also conserving 
momentum and energy 

Discrete Beliaev coupling of excitations was measured 
in @, 0|- Since the spectrum is discrete, such Beliaev 
coupling of low energy excitations is rare, and even at low 
temperatures, damping is usually governed by Landau 
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Damping of bulk excitations 
is due mainly to Beliaev coupling at low temperatures 
11 1 . Such damping was measured [l3|, and found to 
have strong momentum dependence, in agreement with 
However, the dependence of the Beliaev damping on 
the excitation energy was not measured. 

In this letter we study the effects of the damping on 
the excitation line shape in an elongated condensate. 
By comparing the overall response to Bragg excitations, 
with the response of the undamped fraction, the damp- 
ing rates are quantified as a function of both momentum 



and energy of the decaying QP. In addition to the mo- 
mentum dependence of the damping rate, we observe a 
dependence on the excitation ener gy. By modelling the 
condensate as an infinite cylinder [13l|.[14|. we calculate 
the damping rates of the different elementary modes, and 
find good agreement with the measured line shapes. Ac- 
cording to the model, the spatial dependence of the QP 
wave functions plays a key role in the difference between 
the damping rates. The different energy dependence at 
different momentum excitations is a result of quantum 
interference due to many-body effects. 

We create a BEC of iV = 4 x 10^ ^'^Rb atoms in 
the internal state = 2, m^? = 2) in a cylindrically sym- 
metric loffe-Pritchard trap with trapping frequencies of 



27r X 350 Hz and uj^ 



27r X 30 Hz. The Thomas- 



Fermi radii of the condensate are Rtf ~ 3 /xm and 
Ztf ~ 36 ^m, and the chemical potential is « 4 kHz. 
We excite the condensate by shining it with two off res- 
onant laser beams detuned by 0.2 nm from the D2 tran- 
sition for a period of 2ms (Bragg excitation). Then we 
rapidly shut the magnetic trap off, and image the atomic 
cloud after 38 ms time of flight. By varying both the 
frequency detuning lo and the angle a between the two 
beams, we can measure the response of the condensate 
to an excitation with momentum hk = 2hkLSui{a/2) 
in the axial direction z, and with energy hio. Here 
fci = 27r/780nm~^ is the wave number of the lasers. 
Since k^ZTP ^ 1, the wave vector of the excitation k 
is a good quantum number. 

A unique feature of experiments with BEC is the abil- 
ity to observe the damping products in a single image. 
In the insets of Fig. [1] we present time of flight images of 
an excited condensate. In both images roughly one third 
of the condensate atoms were excited to a QP mode with 
momentum 2hkL, but in Fig. [TJd, the damping is much 
larger. This is evident in the plot of the normalized mo- 
mentum distribution in Fig. [T] For the lower energy ex- 
citation (a), the peak in the momentum distribution at 
k = 2kL is more pronounced. For the larger energy ex- 
citation (b), the damping products are dominant. The 
strong reduction of the condensate fraction along with 
the tendency of the damping products toward lower ax- 
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FIG. 1: Beliaev damping of condensates excited with momen- 
tum 2hkL, and different frequencies: 14 kHz (a) and 18 kHz 
(b). Insets - absorption images of excited condensates after 
38 ms time of flight. The graphs show the normalized momen- 
tum distribution p{k) along the axial direction z, which is a 
vertical integral of the absorption images in the insets. The 
free-particle resonance is at 15 kHz and the local density res- 
onance for our condensate is 17.5 kHz. In (a) the undamped 
excitation is larger than that of (b), while the cloud of Beliaev 
damping products is larger in (b). As the Beliaev damping 
extracts atoms from the condensate, the condensate peak is 
greatly suppressed in (b). 



ial momentum are evidence of multiple collisions. 

In order to quantify the damping of the excitations 
at different frequencies we employ a post-selection tech- 
nique [l3l- We measure the response in two ways. First 
we measure the average momentum along the z axis, of 
all atoms in the expanded atomic cloud: (k) /2kL — 
J dk k ■ p{k)/2kL, where p{k) is the momentum distri- 
bution. Since momentum along z is conserved during 
damping, this is a measure of the excitation fraction in- 
cluding damped excitations and will be referred to as 
overall response. The second way to measure response is 
to count the fraction of atoms that remain in the excita- 
tion mode alone: iV2fe^/iVtot, where the occupation N2kL 
and the total number of atoms A^tot are determined by 
a fit to the momentum distribution. This is referred to 
as the undamped population. In the absence of collisions 
the two measuring methods are equivalent. The differ- 
ence between the two measurements is a measure of the 
number of excitations that were damped. 

We repeat this measurement for different values of to. 
In Fig. [5^ we compare the obtained line shape for ex- 
citations with momentum 2hkL. Empty circles are the 
measurements of the overall response, while the filled cir- 
cles are the undamped results of the same images. The 
overall response displays a peak at 16.6 kHz, less than 
predicted by the local density approximation [lij . This 
is due to the relatively large excitation fraction, which 
shifts the resonance towards the free particle value [l^. 
The width of the resonance is due mainly to the inhomo- 
geneous density of the condensate. The response of the 
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FIG. 2: (a) Line shape at fe = 2kL measured in two ways: (o) 
The overall response (fe) /2fc_L, and (•) the undamped pop- 
ulation N2ki^ / Ntot- The dashed and solid lines are gaussian 
fits centered at 16.6 kHz and 15.3 kHz, with widths of 4 kHz 
and 2.6 kHz respectively. The resonance in the line shape of 
the undamped population is shifted down in energy and is 
narrower than that of the overall response. The dotted line 
marks the center of the resonance of the undamped popula- 
tion, (b) The corresponding theoretical curves obtained from 
our model, for the line shape of the overall response (dashed 
line), and that of the undamped population (solid line). The 
radial modes leading to the line shape are marked by bars, 
positioned at their corresponding energy, and with a height 
proportional to their overlap with the condensate. 



undamped population (filled circles) peaks at 15.3 kHz, 
and is 5 times smaller than that of the total population, 
due to Beliaev damping. There is a clear shift down 
in the resonance of the undamped population, implying 
that the damping rate is larger for the more energetic 
excitations. This can be understood intuitively in a local 
density approximation: more energetic excitations are 
in spatial regions of larger density, leading to a higher 
damping rate. 

In Fig. [5^ we present the response of the condensate to 
excitations with smaller momentum - l.lhkL. A reduc- 
tion in the damping as compared to the case of 2hkL is 
evident. This is a result of the quantum interference due 
to the collective nature of the low momentum excitations 
and therefore is present in a homogeneous BEC as well 
[l3 |. A gaussian fit to the overall response is centered 
at 6.48 ± 0.04 kHz. The fit to the undamped line shape 
has a peak at 6.63 ± 0.07 kHz. Contrary to the 2kL case, 
the line shape of the undamped population is not shifted 
down in energy, and the naive linear dependence of the 
damping rate on the local density fails. 

To theoretically account for these effects we include 
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FIG. 3: (a) Line shape at fe = I.I/cl measured both ways: (o) 
The overall response (fc) /l.lfcL, and (•) the undamped popu- 
lation Ni,ikj^ /Ntot- AH circles are averages over 4 data points. 
The dashed and solid lines are gaussian fits giving a central 
frequency of 6.5 kHz for the overall response and 6.6 kHz for 
the undamped population. The dotted line marks the cen- 
ter of the line shape of the undamped population. There is 
no downward shift of the resonance, and no narrowing as are 
present in Fig. [51 (b) The corresponding theoretical curves 
obtained from our model, for the line shape of the overall 
response (dashed line), and the undamped population (solid 
line) . The bars mark the frequency and overlap with the con- 
densate of the relevant radial modes. 



spatial dependence. We exploit the cigar shape of the 
condensate to neglect the z dependence of the ground 
state, thus reducing the problem to ID in the radial di- 
rection [l3| . A similar treatment is performed in [lo| 
to describe Landau damping of quadrupole oscillations 
of an elongated BEC. The ground state V'o(p) of the 
Gross-Pitaevskii Hamiltonian i?GP = -'\/'^/2M + V{p) + 
gN/L IV'ol^ is calculated by imaginary time evolution, for 
a radially dependent potential V{p) — ^Maj^p^ and our 
experimental parameters. Excitations over the conden- 
sate are obtained by solving the Bogoliubov equations: 



huj^u^ (r) 
-hwvVu (r) = 



N 2 

Hgp + IV'ol 

N 2 



Uu + g^iplv^ (1) 

N 2 

Vv + g-^i>l "i^- (2) 



The vifave functions Ui, (r) of Eqs. I1I2I are decoupled to 
(r) — Mn,m,fc(p)e™'^e''^^, and so are Vi, (r). The set of 
quantum numbers v are in this case: k - the momentum 
along z, m - the vorticity around z, and n - the number 
of radial nodes of the wave functions u{p) and v{p). The 
energy of the radial modes hujk,m,n increases with n as 
well as \k\ and Iml. There is a conservation law for k 
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FIG. 4: Damping rates versus the energy of the different radial 
modes for the first few radial modes with k — l.lkL (a), and 
k = 2kL (b). (•) - the calculated damping rate, (o) - the 
damping rate obtained by taking only the terms involving 
u„u„-^u„2 in Eq. [3l 



and m (fc can only decay into q,k — q, and m into mi, m2 
fulfilling mi+m2 = m), but the quantum number n is not 
conserved. We consider the damping of excitations with 
m = as the Bragg pulse carries no angular momentum 
along z. Still, all m values need to be considered as 
damping products. 

The Beliaev coupling is part of the next order 
expansion of the many-body Hamiltonian: Hb — 
X.-u-Mi^b^b^ + H-C. The three QP overlap 



= 2g,/No j dr ( 

+ U„Viy^ Uiy^ + V„Uiy^ V + V^V^^ + V^V^^ V ) (3) 

involves interference of six different quantum pathways. 

Using first order time dependent perturbation theory 
we calculate the damping as a function of time. Since the 
damping is nearly linear, we can extract damping rates 



m,ni ,n2 



,2 sin^ (Acjt/2) 



(4) 



from the different initially excited radial modes v — 
{n,0,k), to all possible pairs of modes vi = (ni,m, g) 
and i>2 — (ji2, — m, k — q), with an energy mismatch of 
All> — bJu — iOy^ — uju2 These rates are plotted in Fig. [4| 
versus the corresponding eigen- frequencies ujn,o,k- 

The radial mode dependence of the damping rate is 
different for the different initial momenta. The damp- 
ing rate from radial modes with momentum fc = 2fci 
increases with n. In this regime, fko 3> /i, excitations 
are nearly single-particle in nature and v (r) 0, for 
both the decaying excitation, and most of the damping 
products. Therefore the damping is mostly governed by 
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the term involving only u,jUp 



This can be seen by 



the fact that the damping rates in Fig.|4)D (filled circles), 
follow the term with only u's (empty circles) [l3|- The 
increase of both rates with radial mode energy is a re- 
sult of the number of available energy-conserving pairs of 
damping products. As the energy of the decaying QP is 
increased, more pairs are available, and the damping in- 
creases. Our calculations show that the n = 4 mode has 
twice as many possible damping products as the mode 
n — 0, in agreement with the damping rate which is ap- 
proximately doubled. 

Damping of modes with momentum l.lfifci is differ- 
ent. First, the energy spacing between modes is com- 
parable to the 2kL case, but the energies are smaller. 
This is manifested in a larger relative increase of avail- 
able modes with radial mode. In fact, the n = 3 mode 
has 5 times more possible damping target modes than 
the n = mode. Alone, this would increase the damping 
from the higher radial modes as happens with the empty 
circles in Fig. |4^. This effect is compensated by interfer- 
ence between the different terms in Eq. [3l Nearly all the 
damping products have ?ia; < /z, and therefore, as the ra- 
dial number increases, the wave functions overlap denser 
regions of the condensate and v (r) becomes more signif- 
icant. Together, the increase due to phase space, and 
the reduction due to quantum interference between the 
terms in Eq. [3] cancel, and the damping is similar from 
all radial modes for k = 1.1 as seen in filled circles of 
Fig.m 

For even smaller momentum excitations, our model 
predicts non-exponential damping. However, the damp- 
ing is so slow, that the infinite cylindrical approximation 
is no longer valid for our trap parameters. 

The line shapes of the overall response, and the un- 
damped population obtained by the model are presented 
below the data in Figs. [2b, Eb- The bars are positioned 
at the eigen-frequencies of the radial modes with the cor- 
responding k (for m = 0). The height of the bars is pro- 
portional to the overlap with the condensate, that deter- 
mines the response to the Bragg pulse 13|. The dashed 



lines are the expected overall response to the Bragg pulse, 
obtained by summation over the responses of the differ- 
ent radial modes including Fourier broadening of each 
mode, as in [l^. The solid line is obtained in a similar 
manner after multiplying each response function by the 
corresponding damping e~^" 0''=*<=ff . We further convolve 
the line shape with a gaussian of width 300 Hz to account 
for residual sloshing in the trap, that changes the effec- 
tive detuning [l^. The suitable effective damping time, 
tef{ = 5 ms, is 3 ms longer than the Bragg pulse. This 
may be due to collisions during the expansion of the con- 
densate and enhancement of the damping due to thermal 
population of the damping products. 

Gaussian fits to the calculated line shapes give res- 
onances at 6.3 kHz for both overall response and un- 
damped population for the I.IA:^ excitations. For the 



2fci excitations, the model gives a line shape centered at 
16.5 kHz for the overall response and at 15.6 kHz for the 
undamped population in good agreement with the mea- 
sured values. The model displays the same narrowing and 
downward shift in the line shape of the undamped popu- 
lation at high momentum, indicating the larger damping 
rate from more energetic excitations. Both model and ex- 
periment also show a uniform, much smaller damping of 
the small momentum excitations. The lowest radial mode 
is more separated from the others, and distinguishable 
from the rest both in the model and in the experimental 
data. 

In conclusion, we use a post selection technique to mea- 
sure the Beliaev damping of QPs in an elongated conden- 
sate. We find a dependence of the damping rate on both 
the excitation energy and the momentum of the decayed 
QP. We measure a shift downward in the excitation en- 
ergy of the undamped excitation for high momentum, 
while for small momentum we measure no such shift. 
Modelling our system as infinite along the axial direc- 
tion, we obtain quantitative agreement between theory 
and experiment. The overall damping is suppressed for 
low momentum excitations as a result of interfering quan- 
tum pathways within each damping channel. 

The inhomogeneous treatment includes a quantization 
of the bath in the radial direction. Since our excitation 
energies are above hujp, we are only indirectly sensitive to 
this quantization. Beliaev damping becomes even more 
intriguing when the radial trapping frequency is increased 
and the system becomes one dimensional. In this regime 
it could be possible to observe effects such as quantum 
Zcno [5| , non exponential damping and even the complete 
inhibition of damping due to the convex excitation line 
shape. 

Studying damping in our system has the advantage 
of imaging which modes of the continuum are excited 
upon damping. This can serve as a spectroscopic tool 
to probe the QP spectrum of the damping products [l9| . 
One can also differentiate between Beliaev damping and 
Landau damping according to the resulting momentum 
distribution of damping products. 

This work was supported by the DIP and Minerva 
foundations. 
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